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Abstract 

Being strictly forbidden in the standard model, experimental detection of the 
lepton flavor violating decays B (B) — > ^~ and 6(6) X fi" would 
constitute an unmistakable indication of new physics. We study these decays 
in supersymmetric models without R-parity and without lepton number. In 
order to derive order of magnitude predictions for the branching ratios, we 
assume a horizontal U{1) symmetry with horizontal charges chosen to explain 
the magnitude of fermion masses and quark mixing angles. We find that the 
branching ratios for decays with a r/i pair in the final state are not particularly 
suppressed with respect to the lepton flavor conserving channels. In general in 
these models B[6 ^ /x+ fi' {X)] < B[6(6) ^ r+ fi' (X)] < B[6 ^ r+ (X)] . 
While in some cases the rates for final states r"*"T~ can be up to one order 
of magnitude larger than the lepton flavor violating channel, due to better 
efficiencies for muon detection and to the absence of standard model contri- 
butions, decays into final states appear to be better suited to reveal this 
kind of new physics. 
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I. INTRODUCTION 



In the standard model (SM) the SU{2) x U{1) gauge symmetry together with Lorentz 
invariance imphes accidental Baryon (B) and Lepton (L) number conservation at the renor- 
malizable level. Due to the larger Lorentz structure, supersymmetric (SUSY) versions of 
the SM allow for renormalizable B and L violating operators involving scalars with non-zero 
B and L charges, that can induce fast proton decay as well as several other unobserved 
processes. Therefore, additional symmetries are required to enforce proton stability and 
to suppress B and L violating transitions. In most SUSY models, invariance under the 
additional parity quantum number R = (_i)3B+l+25 being the spin) is assumed, and 
this enforces B and L conservation at the renormalizable level. However, today it is be- 
lieved that B and L are not likely to be fundamental symmetries of nature, and in fact 
a much larger spectrum of models it is known to be consistent with the data. To render 
phenomenologically viable SUSY extensions of the SM, the first priority is to ensure proton 
stability. In this respect other symmetries can be more effective than R parity, since R par- 
ity still allows for potentially dangerous dimension five B and L violating operators. Some 
interesting alternatives exist which forbid dimension four and five B violating terms, and 
hence are more effective to ensure proton stability [|l|. Since in these models L number can 
be violated by renormalizable operators, they imply a quite different phenomenology from 
R-parity conserving SUSY models 0. 

Two new types of Lagrangian terms characterize this class of models. First we have 
bilinear terms which couple the three lepton doublets superfields , and L^- with the 
up-type Higgs superfield The main effect of these terms is to induce neutrino masses 
via neutrino- neutralino mixing [^]-|^. Requiring that the resulting masses do not exceed 
the experimental limits, implies constraints on the structure of these models [0-^]. In the 
models we will study these constraints are automatically satisfied, thanks to the presence of a 
horizontal symmetry that suppresses all the contributions to neutrino masses. Another effect 
of the bilinear terms is that of mixing fermions in different representations of SU{2). In turn, 
this can generate flavor changing couplings of the Z boson to the leptons. In our theoretical 
framework also these effects are safely suppressed below the experimental sensitivity. For 
these reasons the effects of the bilinear terms do no warrant further elaboration in the present 
context. 

Secondly we have a set of renormalizable interactions in the superpotential which are 
responsible for L and lepton flavor (Lj , z = e , /i , r) violating transitions. In the mass basis, 
these terms read 

Kjk Li Lj -|- X'ijk Li Qj d'j^ , (1-1) 

where Qi and (if denote the quark doublet and down-quark singlet superfields. If are the 
lepton singlets, and Xijk = —Xjik due to the antisymmetry in the SU{2) indices. 

Several of the A and A' couplings are strongly constrained by the existing phenomenology 

0. The best limits are for couplings involving fermions of the first two generations (i, j, k = 

1, 2) while for couplings involving more than a single third generation field the existing limits 
are much weaker. From the theoretical point of view, the values of the A and A' couplings in 
( |L1| ) are not predicted by the model. However, general models that can explain the observed 
fermion mass hierarchy also predict that R-parity violating couplings involving more than a 
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single third generation field are the largest ones. This suggests that this kind of new physics 
can be effectively searched for in rare decays involving b, t and z/^ [0^- In the following, 
we will make this statement more precise, by imposing on the models additional theoretical 
constraints. Following Ref. we embed into the R-parity violating model a particular 
horizontal symmetry that can account for the order of magnitude of the fermion masses 
and CKM angles ||TO|JTT[] . This framework allows us to estimate the size of the relevant L 
violating couplings in (|TTT|). 

A rather complete study of B decays into third generation leptons like B tu^ , B — *■ 
r+r^ , b —>■ Xt^t~ , b — > Xut-Uj. has been recently presented in Ref. [Q. The sensitivity 
of these decay modes to new physics from SUSY models without R-parity was thoroughly 
investigated, and compared to the sensitivity of the corresponding decay modes into /i, i/^. It 
was found that the processes that are most sensitive to these effects are the leptonic decays 
Bd,s T^T~ and the inclusive decay b — >■ X^r+r" . However, from the experimental point of 
view the efficiency for r identification is expected to be rather low at 5-factories. Similarly, 
r-tagging will be a very hard task at future high B statistics experiments at hadron colliders 
|T^. This constitutes a serious drawback for the search of new physics effects in decays with 
final state r's, and it is unlikely that the theoretical enhancement of the decay rates could 
fully compensate for this. 

In this paper we point out that the lepton fiavor violating decays Bd,s{Ba^s) r+yU" and 
b (b) XsT^fi^ together with the corresponding CP conjugate decays, can provide the best 
compromise between the two requirements of large theoretical branching ratios and good 
efficiencies in searching for the experimental signatures. Indeed, for the two body decay 
B T~^fi~ the absence of a with momentum opposite to the /i~ in the B rest frame 
represents a clean signature, rather easy to search for. The first experimental limit on this 
decay B{B r^jJ,^) < 8.3 x 10"^ has been recently established by the CLEO Collaboration 
||r3| . The search for the three body decay b XgT^fi^ appears to be more difficult, because 
of the lack of knowledge of the momentum of the missing /i"*". This is refiected in the present 
experimental situation. While a tight limit on the lepton fiavor violating decay b —>■ X<j/i+e~ 
has been recently established [B(6 — > X^fi^e"^) < 2.2 x 10~^ [14|] to date no experimental 
limit exists on decays into XgT'^e" or XgT^fi^ final states. 

At hadron colliders, already the study of the two body decay B r^fi^ will be a difficult 
task. This is because in this case the momentum of the decaying B is not known, hence the 
presence of large backgrounds, as for example from the semileptonic decays B D'^*^ iiu ^ 
will render quite challenging the search for this rare decay. 

From the theoretical point of view, the detection of lepton fiavor violating decays would 
represent a striking evidence of physics beyond the SM. The absence of SM contributions, 
and in the case of the decay b XsT^ ^~ the absence of long distance effects which are 



difficult to estimate in a reliable way [|I5|, render these decays well suited to reveal in a 
clean way new physics effects. In summary, because of the large theoretical enhancement 
of the branching ratios with respect to yU'^/i" final states, and since in any case a muon is 
experimentally much easier to identify than a tau, we believe that these processes will allow 
to search for signals of SUSY models without R parity with a better sensitivity than lepton 
flavor conserving decays into fi~^fi~ or t^t~ . 

In the next section we will outline the main features of SUSY models without R-parity 
embedded in models with horizontal symmetries. In Sec. Ill we will present the relevant 
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expressions for the effective new physics coefficients which appear in these models, and for 
the decay rates. Finally, in Sec. IV we will discuss our results and present our conclusions. 



II. R-PARITY VIOLATION IN THE FRAMEWORK OF HORIZONTAL 

SYMMETRIES 



In order to evaluate the effects of the R-parity violating interactions in ( |1.1| ), we need 
to estimate quantitatively the coefficients A and A'. We work in the framework of the 
supersymmetric models with horizontal symmetries that have been thoroughly investigated 
in Refs. ||To| , p!T|| . We assign to each supermultiplet ip a charge H{ip) of an Abelian horizontal 



group Ti. = U{1)h which is explicitly broken by a small parameter e with charge H{e) = —1 . 
This gives rise to a set of selection rules for the effective couplings appearing in the low 



energy Lagrangian |]10[. Assuming that each of the lepton, quark and Higgs superfields 
carries a positive or zero charge, the selection rule relevant for the present discussion is 
that the effective coupling gate for a general trilinear superpotential term ipai^bipc is of order 
dabc ~ £H{ipa)+H(ip,,)+H{tpc) _ Therefore, the leptons and down-type quarks Yukawa couplings 
are respectively of order YA ~ ^H(^,)+HiL,)+Hii^,) and Yf- ~ . Most of 

the L-violating couplings in ( p. . 1| ) are further suppressed with respect to the corresponding 
Yukawa couplings. They can be estimated as 

1 /2 

A,., ~ yA ^ f , (2.1) 

\ COS l) j 

and 

1 /2 

J J \^ cos^ p y 

where Gp is the Fermi constant, and tan {3 = with the up-type Higgs doublet. 

From Eqs. (|2.1| ) and (|2.2| ) it is apparent that in our framework the couplings A and A' involving 
fermions of the third generation are respectively enhanced by m,- and rrib ■ 
In order to give a numerical estimate of the couplings, we need a set of H charges and a 
value for e . The magnitude of the Ti. breaking parameter is generally taken to be the value 
of the Cabibbo angle, e ~ 0.22 , while the quark, lepton and Higgs charges are chosen to 
reproduce the values of the fermion masses and CKM mixing angles. Besides reproducing 



the measured values, the model has some predictivity in the quark sector ||10[, it yields 
estimates for ratios of neutrino masses pJ] , p!6[] and, most important in the present context, 
it ensures that the L-violating couplings in (|L1|) , (p.l|) and ( p.2|) are safely suppressed below 
the present experimental limits [Q. The following ?i-charge assignments fit the order of 



magnitude of all the quark masses and CKM mixing angles ||T0|| : 



Qi Q2 Q3 dl dl dl ill ul ul $d 
(3)(2)(0) (3) (2) (2) (3)(1)(0) (0) (0) . (2.3) 

Following Ref. 0, we use for the leptons two different sets of charges which define two 
different models, and for each model we chose a different value of the squark masses rriq : 
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Li L2 U k k k m;(GeV) (GeV) 

Model I : (4) (2) (0) (4) (3) (3) 100 170 

Model II : (3) (0) (0) (5) (5) (3) 100 350 . (2.4) 

Model I tends to enhance new physics effects induced by operators arising from squark 
exchange, while in model II the effects of new scalar operators induced by slepton exchange 



tend to dominate. The choices ( p. 3D and (|2.4|) for the horizontal charges are not unique. 



Since the Yukawa interactions are invariant under a set of U{1) symmetries such as B, L 
and hypercharge, it is always possible to shift the if-charges of any amount proportional to 
one of the corresponding f/(l) quantum numbers without affecting the predictions for the 
masses and mixing angles. In particular, the shift (proportional to L) H{Li) — > H{Li) +n, 
H{lf) H(lf) — n and H{ip) — > H{ip) for all the other fields has the effect of suppressing 
(for n > ) all the L violating couplings in (pTTp , (|2.1| ) and (|2.2|) by a factor of . It turns 
out that already for n = 1 the suppression is large enough so that all the lepton flavor 
violating decays will be unobservable at most of the future -B-physics experiments. We also 
notice that model II can be derived from model I by means of shifts proportional to lepton 
flavor numbers with Uf, = —1 , = —2 , = . This has the effect of enhancing some of 
the A couplings without affecting the charged lepton masses. Of course, the predictions for 
the neutrino mixing angles will be different in the two models. 



III. COEFFICIENTS AND RATES FOR THE DECAYS 

In the models under investigation, the lepton flavor violating decay Bq — > r+/x~ 
(with q = d,s) and b — > XsT^fi^ are induced by the effective Lagrangian 

- Cff = {qhhR) {UrTl) + C^s {qBbh) {^lTr) 

+ Cy {qRi^hR) {fiLl^^rL) + h.c. . (3.1) 

The first two operators arise from sneutrino exchange diagrams, while the last one corre- 
sponds to squark exchange diagrams after Fierz rearrangement. The coefficients appearing 
in (|3.1|) read 

where the index values i = 3 in C^g and i = 2 in are excluded because of the antisym- 
metry in the first two indices of the A couplings. 

For the decays Bq T~^fi~ and b X^r'^fi' the effective Lagrangian reads 

- ^ts = Cts ibLqR) {UrTl) + Cts (bRqi) {UlTr) 

+ C+ (hR^^qR) ifLLl^^TL) + h.c. , (3.3) 



where 



5 



The expressions for the various branching ratios are presented in the next two subsections. 
In order to simphfy the formulae we have neglected the muon mass (however, 7^ has 
been kept in the numerical analysis). 

A. The decays B —>■ r+/i and B —>■ r^fi^ 
The amplitude for the decay Bg r^fi^ can be written as 



A =ifBmB- 



(fir) 



+ 



rriB 
nib 



i^2S + ^i; 



m 



TUB 



(/^75i-) 



(3.5) 



where we have used the PCAC (partial conserved axial current) relations (0| wf^'-f^b \B) 
ifsPB and (0| U75 6 \B) ~ —ifs ^'s/^b ■ This yields the branching ratio 



BiBr, 



T H 



e2 



fsjB 



ml 



64 TT 



m 



B , 



rriB 



rrir 



^2S 

rrib TUB 



a: 



+ 



m 



B 



mf 



(3.6) 



Eq. ( p.6|) also accounts for the CP conjugate decay Bg yU+ , therefore experimental 

searches for both the decay channels B r+yU^ and B T^fi~^ will yield informations 
on the same set of operators. The decay mode B —>■ r^fi^ is controlled by the coefficients 



C{s, C^s and in 



The amphtude is given by 



-AziCrs,C^s,C^) with A- 



defined as in ( p.5|) . Therefore the branching ratios B(B ji ) and B(i? r /i"'') are 

again given by ( p. 61) with the substitution {Cf^, , Cy} {Cts^ ^ts ' ^v}- 



B. The decays h Xt^ij, and b Xt^/j, 



For the double differential distribution for the inclusive decay b{p) — > s{p') T~^{k') ^ (k) 
with respect to the invariants x = p' ■ k/ml and y = p ■ k/ml , we find 



mr 



dxdy IGir^ 



4|a- 



1 — mi — mz 



- |2 
15 1 



x-y 



y y 



1 — mz. + m„ 



(x-y) - 2Re{C2sCy*)mrX 



where rhf = mj/mb, and x and y range between the following limits 

1 — (rhs + m^)^ 



0<y< 



X- < X < x+ , 



with 



x± 



2s 



m. 



+ (1-S)( 



- 2 I \ 1/2/ - - 2 - 2- 

s — m^ ± A ' [s, m^, m^ 



(3.7) 



(3.8) 



(3.9) 



where s = {p — p') /m^ and A(x, y, z) = x + y + z — 2xy — 2xz — 2yz. 
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The corresponding expressions for the distribution d?T^ /dxdy describing the decay h — > 
s r+ /i~ can be obtained from ( p.7|) by interchanging p ^ —p' which yields x ^ —y and 
^ 1 in the terms inside square brackets, and by substituting {Cf^, Cg^jCy} with 

We now introduce the forward-backward asymmetries Apg of the two distributions . 
The asymmetries are defined with respect to the angular variable ce = cos 6, where 6 is the 
angle between the /i~ momentum k and the s momentum p' in the B rest frame: 



dT^y)/dy 



dco : — : / dCf 



dydce J-i dydc^ 



(3.10) 



As we will see, the dependence of the asymmetries with respect to the normalized muon 
energy y = E/j^/rrib can provide important informations on the underlying new physics, 
which are complementary to the measurements of the branching ratios. The average values 
of the asymmetries are computed as 

(^Ib) = ^ / dy^^^ A$M , (3.11) 

where = / dy '^^^^^^ . Finally, we also study the total /i^ asymmetry ApB of the distribu- 
tion d"^ [F~(?/, cg) + F+(?/, co)]/dydcs which turns out to be a useful quantity when untagged 
B samples are used for the measurements. 



IV. DISCUSSION 

In this section we discuss the numerical predictions for the branching ratios for the 
decays B^^g (-B^,.)) — T~^fJ'~, b (b) — > XsT~^fj,~ and for the /i~ forward-backward asymmetries 
measurable in the decays into three body final states. In our estimates, we have used the 
following set of values for the relevant SM parameters: = 4.8 GeV, = 200 MeV, = 
106 MeV, TTir = 1.777 GeV, /s, = 200 MeV, fe, = 230 MeV, = 5.3GeV, tb = 1.6ps 
and B{B Xclv) = 10.4% , while the magnitude of the various new physics coefficients is 
determined by the sets of if-charges and SUSY masses listed in ( ^l3|) and ( |2.4|) . 

Our results are collected in Table I and in Figs. 1-3. Table I lists the numerical results 
for B decays involving the channels b —>■ fi~ and b —>■ fi~ , which are respectively controlled 
by the two different sets of coefficients {Cig2sv} {^'1525^}- These coefficients are 
evaluated in the two different models defined by the charges in ( |2.3| ) and ( p.4| ). The entries 
in the first column in Table I refer to model I. We recall that in this case the choice of the 
leptonic horizontal charges tends to enhance the effects of squark exchange. The entries 
in the second column refer to model II. Here squark exchange diagrams are suppressed by 
a different choice of the horizontal charges and by the relatively large value of the squark 
masses, so the effects of slepton exchange tend to dominate. We stress that the aim of the 
numerical predictions given in the first two columns in Table I is that of suggesting the 
level of precision that future -B-physics experiments will have to reach in order to detect, or 
to effectively constrain, new physics from SUSY without R-parity. In the last column we 
list the existing experimental limits. It is apparent that most of the decay modes we have 
studied are presently still unconstrained. 
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TABLE I. Predictions for the branching ratios for the lepton flavor violating B and B de- 
cays into ^~ and for the forward-backward /x~ asymmetries, in the R-parity violating models 
discussed in the text. Model I is defined by the lepton horizontal charges H{L) = (4,2,0), 
= (4,3,3) and by the SUSY masses mj = 100 GeV, rrig = 170 GeV. Model II corre- 
sponds to the horizontal charges H{L) = (3,0,0), H(l'^) = (5,5,3) and to in^ = 100 GeV 
and rriq = 350 GeV. In both models the value of the horizontal symmetry breaking parameter is 
e = 0.22 . The existing experimental limits are given in the third column. 



Process 


Model 1 


Model 2 


90% c.l. hmit [Ref.] 


B {Bs ^ T+fi') 


8.3 X 10-^ 


7.9 X 10-^ 






B (Bd ^ r+/i-) 


3.0 X 10-1° 


2.9 X 10~^ 


8.3 X 10"^ 


m 


B {B, ^ T+fi-) 


5.0 X 10"^ 


2.7 X lO^'' 






B (Bd ^ T+fi-) 


1.8 X lO^s 


1.0 X 10-5 


8.3 X 10-^ 


131 


B {b^XsT+fl-} 


1.9 X lO"'^ 


6.4 X 10-5 






B (b^XsT+i^-) 


1.6 X lO"'^ 


4.1 X 10~6 






(^?B (^)> 


-0.02 


-0.24 








-0.75 


-0.71 






{Afb {b,b)) 


-0.36 


-0.27 







The first four lines in Table I collect tlie results for tlie two body leptonic decays, wliile tlie 
results for the three body final states are given in the following two lines. Next we present the 
results for the average values of the /i~ forward-backward asymmetries Ap-^ corresponding 
respectively to 6 — > /i~ and b ^ fi~ decays. The average value of the untagged fi~ asymmetry 
ApB is given in the last line. 

To put in evidence the advantage of studying the lepton flavor violating decay modes 
with respect to the lepton flavor conserving decays b r^T^{X) and b fi^fi~{X) , we 
list in Table II (taken from Ref. 0) the numerical predictions for these decays as derived in 
the SM, in model I and in model II. 

Our results for the two body decays are as follows: the rates for B decays are more 
than two orders of magnitude larger than the corresponding B decays. This is mainly 
due to the fact that the dominant contribution to the two body decays of B mesons is 
given by the coefficient , while for the B , is suppressed by the if-charge difference 
H{d'^) — H{Q^) = 2, so that it gives negligible contributions to the decay rates. This is 
interesting, because - modulo generation independent shifts proportional to baryon number 
or hypercharge - the set of charges in (|2.3| ) is unique for fitting the quark masses and mixing 
angles. Therefore, the hierarchy in the pattern of the decay rates shown in Table I can 
be taken as a general quahtative prediction of SUSY models without R-parity embedded in 
models with an U{1) horizontal symmetry. Confronting the figures of the leading Bg — > r+/i~ 
decay rates with the corresponding figures in Table II, we see that while in model I the rates 
are much smaller than the rates for Bq t~^t~ , in model II there is no similar suppression. 
In this case the larger phase space available for the lepton fiavor violating decays yields 
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TABLE II. Predictions taken from Ref. Q for the lepton flavor conserving B decays into r+r" 
and in the standard model and in the presence of new physics. Model I and model II 

coincide with the two models of Table I, and are discussed in the text. The existing experimental 
limits are listed in the last column. 



Process 


Standard Model 


Model 1 


Model 2 


Limit 




B {Bs T+T-) 
B (Brf ^ T+T-) 


9.1 X lO"'^ 
4.3 X 10-^ 


5.7 X 10-6 
1.9 X lO-'^ 


1.8 X lO-'' 
6.3 X 10-6 


5.0 X 10-2 ^ 
1.5 X 10-2 ^ 


i 
i 


B (Bs ^ 
B {Bd ^ 


4.3 X 10^9 
2.1 X 10-1° 


7.9 X 10-^ 
2.9 X 10-^ 


7.2 X 10-^ 
2.7 X 10-^ 


2.6 X 10-6 ^ 
8.6 X 10-^ ^ 


171 
17] 


B (6 ^ Xsr+r")nocut 
B (6^X,r+T-),>o.6 


4.9 X 10-6 
1.5 X lO-'^ 


7.3 X 10-6 
2.2 X 10-6 


7.9 X 10-6 
2.7 X 10-6 


5.0 X 10-2 ^ 


i 


B (6 ^ 

B (6 ^ Xs fi+fi')i^oA 


3.1 X 10-"^ 
4.3 X 10-6 


3.1 X 10-"^ 
4.4 X 10-6 


3.4 X 10-"^ 
8.4 X 10-6 


5.8 X 10-^ 1 


14 1 



Limit estimated from the non-observation of large missing energy events at LEP. 



^ 95% C.I. " 90% C.I. 



a slight enhancement of the rates. Confronting with the experimental limit on this decay 
mode given in the third column in Table I, we see that an improvement of two orders of 
magnitude is needed in order to test the model. Such an improvement can be within the 
reach of forthcoming 5-factory experiments. 

In contrast to the two body decays, the decay rates for b and b into three body final 
states predicted in model I are comparable in size. This is because the leading contribution 
to both decays now comes from vector operators, and confronting ( p.2|) with (|3.4|) it is easy 
to check that Cy = Cy. Model II tends to enhance the effect of scalar operators, and 
this again results in a relative enhancement of the b with respect to the b decay rates. A 
confrontation with the figures given in Table II shows that in model II the lepton flavor 
violating b decay can be up to one order of magnitude larger than the rate for b — >■ XsT~t~^. 
Again, this is mainly due to phase space effects. Table II also shows that without kinematic 
cuts, the SM rates for three body final states are comparable with the rates predicted by 
the new physics models. Then, in order to single out the new physics short distance effects 
in the lepton flavor conserving channels, it is necessary to impose suitable cuts. Table II 
shows the effects of two different cuts on the t~^t~ and fi~^fi~ invariant masses. Clearly, one 
of the advantages of studying lepton flavor violating channels is the complete absence of this 
background to the new physics effects. 

In spite of some large enhancements from the new physics, the figures in Table I make ap- 
parent that even in the most favorable cases, the predicted branching ratios remain rather 
small. Therefore, in order to measure the corresponding rates, or to put significant con- 
straints on the models, a large statistics and a good experimental efficiency are required. It 
is worth noticing that while a separate measurement of the relevant combinations of the new 
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physics coefficients {Cj^25y} ^^^1 {C'^gsy} requires the identification of the flavor of the 
decaying B , with a corresponding loss in the experimental efficiency, in order to establish 
limits on the new physics coefficients it is sufficient to search for r/x pairs production in the 
decays of untagged B samples. This procedure can ensure a large gain in statistics and will 
yield the strongest bounds. If t/j, pairs are detected in the decay of a sample of untagged 
S's, then a measurement of the fi forward-backward asymmetry could provide the additional 
information needed to identify the flavor of the initial state. 

Figs. 1 and 2 depict respectively the forward-backward asymmetries for b and b as 
a function of the normalized muon energy in the B rest frame y — E^/ rrih . The solid lines 
correspond to model I, while the dashed fines refer to model II. We see that while the shape 
of the asymmetries is quite similar for the two models, there are large differences in the 
H~ angular distributions for 6 or 6 initial states. Namely, the asymmetry for decaying b is 
negative in the whole energy range, yielding the large negative averages listed in Table I. In 
contrast, for decaying 6's the asymmetry changes sign. This induces large cancellations in 
the energy averages as is apparent from the figures in the next-to-last fine in Table I. Fig. 3 
depicts the asymmetry for an untagged sample of an equal number of b and b initial states. 
Since in model II (dashed line) the decay rate for b — » XsT^fi^ dominates over the rate rate 
for b , the asymmetry for the untagged sample reproduces quite closely the energy dependence 
of the asymmetry in Fig. 1 (namely there is a change in the sign) . In model I the two decay 
modes have comparable branching ratios, and accordingly it is not possible to identify in a 
reliable way the flavor of the initial state just from an inspection of the untagged asymmetry 
(solid line in Fig. 3). From these results it is clear that a measurement of a change of sign 
in the energy dependence of the /i^ asymmetry in an untagged B sample would signal that 
most of the r+/i" events originates from decays of 6's, while the measurement of a negative 
asymmetry over the whole energy range would suggest that the contribution from decaying 
6's is at least comparable in size. 

To check to what extent this remarkable feature of the asymmetries depends on our 
particular models, we have arbitrarily varied the values of the scalars and vector coefficients 
in the two models. Our results indicate that in the limit of very heavy squarks rriq 3> 
350 GeV, which implies Cy <S Cg , the difference in the energy dependence of the two 
asymmetries tends to be washed out: namely also for initial b we find a sign inversion. In 
the opposite limit Cf <^ Cy , corresponding to very heavy sleptons, both the asymmetries 
become negative over the whole energy range, even if there are large differences in the 
shapes. However, if the difference between the squarks and sleptons masses is kept within 
a few hundreds GeV, the qualitative features of the energy dependence of the asymmetries 
are maintained, rendering possible in principle the identification of the flavor of the initial 
state. 

In conclusion, in this paper we have investigated b and b lepton flavor violating decays 
involving a r/x pair in the flnal state. These decays are expected to occur in SUSY models 
without R-parity. Assuming a L'^(l) horizontal symmetry with fermion charges chosen to 
explain a known set of parameters ( the quark masses, the CKM mixing angles, and the 
lepton masses) allowed us to derive order of magnitude predictions for the various branching 
ratios. A straightforward prediction of this theoretical framework is that new physics effects 
are stronger in decays like b — > r+r^(X) when several third generation fermions are involved. 
However, our results indicate that in general decays involving a t/j, pair in the final state 
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are not particularly suppressed with respect to the decays involving a pair of r's. On the 
other hand, due to the presence of a single muon in the decay products and to the absence 
of any SM contribution (and in particular of the backgrounds from long distance effects) the 
lepton flavor violating channels are experimentally much easier to be searched for. Therefore, 
the decays B{B) — > r"^//" and b(b) — > XsT'^ji" , together with the CP conjugated decays, 
appear to be better suited than the lepton flavor conserving decays to search for signals of 
R-parity violation. 
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FIG. 1. Predictions for the /x~ forward-backward asymmetry Ap^{y) in the decay b ^ Xg T~^fi~ 
as a function of the normahzed B rest frame muon energy y = Efj^/mf,. The sohd (dashed) 
hne correspond to model I (II) discussed in the text. Model I is defined by the lepton hori- 
zontal charges H{L) = (4,2,0), i7(r) = (4,3,3) and by the SUSY masses m^- = 100 GeV, 
rriq = 170 GeV. Model II corresponds to the horizontal charges H(L) = (3,0,0), H(l'^) = (5,5,3) 
and to nil ~ ~ GeV. In both models the value of the horizontal symmetry 

breaking parameter is e = 0.22 . 
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FIG. 2. Predictions for the /x^ forward-backward asymmetry Ap^{y) in the decay b ^ Xg t~^ij,~ 
as a function of the normalized B rest frame muon energy y = E^/nii,. The solid (dashed) line 
correspond to model I (II) discussed in the text. The new physics model parameters are as in 
Fig. 1. 
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FIG. 3. Predictions for the forward-backward asymmetry ApB{y) in the decay of an un- 
tagged b sample into Xg T'^fi~ as a function of the normahzed B rest frame muon energy y = Efj^/mi,. 
The sohd (dashed) hne correspond to model I (II) discussed in the text. The new physics model 
parameters are as in Fig. 1. 
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